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Micromechanics  of  An  Extrusion  in 


* 

High-Cycle  Fatigue 

by 

T.H.  Lin  ,  S.R.  Lin  ,  and  X.Q.  Wu 


Abstract 

^  A  most  favorably  oriented  crystal  located  at  a  free  surface  of  a  f.c.c.  polycrystal  under 
cyclic  tension  and  compression  of  high-cycle  loading  is  considered.  An  extrusion  in  this  cry¬ 
stal  is  shown  to  be  caused  by  a  positive  slip  in  one  thin  slice  ”P"  and  a  negative  slip  in  a 

closely  located  slice  "Q".  An  initial  tensile  strain  e^  in  the  thin  slice  "R"  sandwiched 

\ 

between  P  and  Q  causes  an  initial  compressive  stress  z^a  in  R  and  a  positive  initial  shear 
stress  in  P  and  a  negative  one  in  Q.  (The  repetition  of  the  subscript  in  Greek  letters  in  this 
paper  as  the  above  does  not  denote  summations).  Slices  P,  Q  and  R,  slip  direction  a  and  nor¬ 
mal  to  the  slip  plane  P  all  make  45°  with  the  free  surface.  The  elongation  in  R  induced  by 
this  initial  strain  is  called  the  "static  extrusion".  The  difference  in  resolved  shear  stresses  in  P 
and  Q  causes  the  build-up  of  plastic  shear  strain  in  P  and  Q,  hence  the  extrusion  growth.  As 
the  extrusion  grows,  the  initial  compression  in  R  decreases  resulting  in  a  decrease  in  the 
extrusion  growth  rate.  This  decrease  of  compression  in  R  tends  to  activate  a  second  slip  sys¬ 
tem  to  slide.  The  plastic  strain  due  to  slip  in  this  second  slip  system  has  a  tensor  component 
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eaa ,  which  has  the  same  effect  as  the  initial  strain  e‘aa  in  causing  the  positive  and  negative 
resolved  shear  stresses  in  P  and  Q,  and  hence  the  additional  extrusion  growth.  In  the  present 
study  with  the  consideration  of  the  secondary  slip  in  R,  the  extrusion  is  shown  to  grow  far 
beyond  the  static  extrusion.  A  similar  conclusion  can  be  made  for  the  growth  of  an  intrusion. 

Introduction 

Forsyth  (1953)  made  an  important  discovery  of  extrusion  and  intrusion  in  fatigue 
bands  in  aluminum  alloys.  Thompson,  Wadsworth  and  Louat  (1956)  and  Hall  (1958)  detected 
extrusions  in  both  copper  and  aluminum.  Extrusions  were  observed  by  Meke  and  Blochwitz 
(1980)  and  Mughrabi  (1980)  in  their  studies  of  persistent  slip  bands.  Extrusions  and  intru¬ 
sions  in  fatigue  specimens  were  also  observed  by  a  number  of  other  investigators. 

Following  the  clue,  provided  by  the  observation  on  extrusions  and  intrusions,  a  number 
of  theories  of  fatigue  crack  initiation  have  been  proposed  by  different  distinguished  investiga¬ 
tors:  Mott  (1958),  Cottrell  and  Hull  (1957),  Thompson  (1959),  McEviley  and  Machlin  (1959), 
Wood  (1956)  and  others.  One  theory  considered  a  column  of  metal  containing  a  single  screw 
dislocation  intersecting  a  free  surface.  Assuming  that  this  dislocation  travels  a  complete  cir¬ 
cuit,  the  volume  contained  in  the  circuit  is  translated  parallel  to  the  dislocation.  This  causes 
the  metal  to  extrude.  However,  this  mechanism  does  not  explain  why  the  dislocation  under 
cyclic  stressing  does  not  oscillate  back  and  forth  along  the  same  path  rather  than  transversing 
a  closed  circuit.  Clearly,  some  form  of  gating  mechanism  is  required  to  convert  the  back  and 
forth  oscillation  of  the  dislocation  into  a  uni-directional  circuit.  This  gating  mechanism  is 
also  needed  in  a  number  of  other  proposed  mechanisms  a«  discussed  by  Kennedy  (1963). 
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It  is  well  known  that  single  crystal  tests  show  that,  under  stress,  slip  occurs  along  cer¬ 
tain  directions  on  certain  planes;  this  slip  depends  on  the  resolved  shear  stress  and  not  on  the 
normal  stress  on  the  sliding  plane.  This  dependency  of  slip  on  the  resolved  shear  stress, 
known  as  Schmid’s  law,  has  been  shown  by  Parker  (1961)  to  hold  also  under  cyclic  loadings. 
Hence  to  calculate  the  slip  distribution  in  a  fatigue  band,  we  need  to  calculate  the  resolved 
shear  stress  distribution  in  the  band. 

Referring  to  a  set  of  rectangular  coordinates,  we  consider  the  strain  to  consist  of  the 
elastic  part  ctJ  and  the  inelastic  part  e^: 

eij  =  eij  +  eij 

(1) 

The  inelastic  strain  may  be  a  combination  of  thermal  strain,  ej,  creep  strain  e^,  plastic  strain 
and  initial  strain  e,j.  Consider  the  elastic  constants  to  be  isotropic.  The  stress  is  related  to 
the  elastic  strain  as 

^ij  —  ^ij  ^kk  +  e,j 

(2) 

where  X  and  G  are  Lame’s  constants,  5;j  is  the  Kronecker  delta  and  the  repetition  of  subscript 
denotes  summation. 


The  condition  of  static  equilibrium  is  written  as 
'■ F,  ~  0  ;  ^ijvj  =  Tiv 

(3) 

where  the  subscript  after  comma  denotes  differentiation  with  respect  to  the  coordinate  vaii- 
able;  F;,  is  the  i-th  component  of  the  bodv  force  and  TjV  is  the  i-th  component  of  the  «urfae? 
traction  on  the  surface  with  outward  unit  normal  v. 
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Substituting  Eq.  (2)  into  Eq.  (3)  yields 


+  2  Geyj  +  F,  -  (Xe^;  +  2G  e^)  =  0 
(Xe^  Sjj  +  2  Gey)  v,-  =  T,v  +  (Xe^  Sy  +  2  Ge")  Vj 


(4) 

(5) 


It  is  seen  that  -(Xe^y  +•  2Ge”j)  is  equivalent  to  F,  and  (Xe^  5y  +  2Gey)Vj  is  equivalent  to  T,v 
in  causing  the  strain  field  e^;  hence  they  are  called  equivalent  body  and  surface  forces  and  are 
denoted  by  Fj  and  T,v  respectively.  The  stress  field  is  then 


^ij  ~  Sjj  X  (e^  e^j)  +  2G(ey  ey ) 


(6) 


Hence,  the  strain  distribution  in  a  body  with  inelastic  strains  under  an  external  load  is  the 
same  as  that  in  a  purely  elastic  body  without  inelastic  strains  but  with  additional  equivalent 
body  and  surface  forces.  This  analogy  reduces  to  Duhammel’s  analogy  of  thermal  stress  (Lin, 
1968)  when  the  inelastic  strain  is  thermal  strain  and  gives  the  same  results  as  those  given  by 
Eshelby’s  ingenious  process  of  cutting,  relaxing,  restoring,  welding,  and  relieving  in  his 
famous  paper  on  ellipsoidal  inclusions  (Eshelby,  1959).  This  method  of  equivalent  force  is 
here  used  to  calculate  the  stress  fields  in  fatigue  bands. 


A  Polycrystal  Model 

Fatigue  cracks  generally  initiate  at  a  free  surface.  To  relieve  the  same  amount  of 
resolved  shear  stress  in  a  thin  slice,  a  greater  amount  of  slip  i.e.,  plastic  strain,  is  required 
near  the  free  surface  than  at  the  interior  of  metals  (Lin  and  Ito,  1967).  For  the  present  analyt¬ 
ical  study,  the  thin  slices  subjected  to  alternate  forward  and  ^ve-ccd  loadings  are  taken  to  he 
in  the  most  favorably-oriented  crystal  at  the  free  surface  of  a  polycrystal.  The  slip  plane  and 
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the  slip  direction  of  the  most  favorable  slip  system  of  this  crystal  form  45°  with  the  specimen 
axis.  The  polycrystal  is  subject  to  a  cyclic  tension  and  compression  of  low  amplitude,  so  the 
plastic  deformation  essentially  occurs  only  in  this  most  favorably  oriented  crystal. 

Defects  exist  in  all  metals.  The  initial  resolved  shear  stress  field  caused  by  initial 
defects  is  denoted  by  x1.  During  loading,  when  the  critical  shear  stress  is  exceeded  in  some 
region,  slip  occurs.  After  unloading,  this  slip  remains  and  induces  a  residual  resolved  shear 
stress  field  xr.  The  resolved  shear  stress  caused  by  loading  is  denoted  by  xa.  Hence,  the  total 
resolved  shear  stress  is  the  sum  of  these  three  stresses:  x1  +  xr  +  xa. 

For  an  extrusion  to  initiate,  positive  shear  has  to  occur  in  a  thin  slice  P  and  negative 
shear  in  a  closely  located  slice  Q,  as  shown  in  Figure  1.  The  initial  stress  field  x1  favorable 
for  this  slip  sequence  is  one  having  positive  shear  stress  in  P  and  negative  in  Q.  Lin  and  Lin 
(1983)  have  shown  that  such  an  initial  stress  field  can  be  provided  by  an  initial  strain  in  the 
thin  slice  R,  sandwiched  between  P  and  Q.  Referring  to  Fig.  1,  an  initial  strain  field  e^a 
varying  linearly  from  zero  at  the  free  surface  to  a  maximum  value  at  the  interior  boundary  of 
the  crystal,  causes  a  constant  gradient  of  strain  e„a  along  the  a  direction  and  a  uniform 
equivalent  body  force  Fa  in  R  yielding  a  uniform  initial  positive  resolved  shear  stress  xl  in  P 
and  a  negative  one  in  Q.  Lin  and  Ito  (1969)  suggested  that  this  positive  e^,  can  be  produced 
by  a  row  of  interstitial  dislocation  dipoles,  and  a  negative  e^a  by  vacancy  dipoles.  Recently, 
Antoncpoulus,  Brown  and  Winter  (1976)  and  Mughrabi,  Wang,  Differt  and  Essmann  (1983) 
have  shown  that  the  ladder  structure  in  a  persistent  slip  band  (PSB)  can  be  represented  by  an 
array  of  dislocation  dipoles  causing  initial  resolved  shear  stresses  at  the  interface  between  the 
PSB  and  the  matrix. 
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Fig.  1  Most  Favorably  Oriented 


Crystal  at  Free  Surface 
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Gating  Mechanism  Provided  by  Micro  Stress  Fields: 


Consider  x1  to  be  positive  in  P  and  negative  in  Q  as  that  caused  by  an  aiTay  of  intersti¬ 
tial  dislocation  dipoles  discussed  earlier.  Referring  to  Fig.  1,  a  tensile  loading  causes  a  posi¬ 
tive  resolved  shear  stress  xa  in  this  crystal.  In  P  the  shear  stress  will  be  the  sum  of  x,  and  xa 
and  will  be  the  first  to  reach  the  critical  stress  xc  and  slide  first.  The  plastic  strain  and  the 
equivalent  forces  are  taken  to  be  constant  along  the  x3  direction,  5xa3/dx3  =  0.  From  the 
equilibrium  condition  Eq.  (3),  with  no  body  force  we  have  3xaa/dot  +  dx^/Bp  =  0,  where  p 
is  normal  to  the  slip  plane.  Since  3xact/3 a  is  finite,  3xaj3/8p  is  also  finite  and  the  change  of 
xap  across  the  short  distance  between  P  and  Q  is  very  small.  Therefore,  the  slip  in  P  relieves 
not  only  the  positive  shear  stress  in  P,  but  also  in  its  neighboring  region  including  Q  as  shown 
by  Lin  and  Ito  (1969)  and  Lin  (1977).  This  keeps  the  positive  shear  stress  in  the  neighboring 
region  from  reaching  that  of  P,  hence  only  P  slides  during  the  forward  loading.  This  slip 
increases  the  negative  resolved  shear  stress  in  Q  to  cause  Q  to  slide  more  readily  in  the 
reverse  loading.  When  the  negative  resolved  shear  stress  in  Q  reaches  xc  in  the  reversed  load¬ 
ing,  slip  occurs  in  Q  and  a  new  residual  stress  field  is  produced.  Similarly,  the  slip  in  Q 
causes  the  relief  of  negative  resolved  shear  stress  not  only  in  Q  but  also  in  its  neighboring 
region  including  P,  thus  making  P  more  readily  to  slide  in  the  next  forward  loading.  This 
process  is  repeated  for  every  cycle,  providing  a  natural  gating  mechanism  for  the  monotonic 
buildup  of  local  slip  and  plastic  strain  in  P  and  Q.  The  plastic  shear  strains  eap  in  P  and  Q 
tend  to  push  out  the  slice  R  and  start  an  extrusion.  Interchanging  the  signs  of  initial  stresses 
in  P  and  Q  will  initiate  an  intrusion  instead  of  an  extrusion.  This  theory  is  extensively  sup¬ 
ported  by  metallurgical  observations  (Lin,  1977).  Some  of  them  are  given  below. 
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Some  Related  Experimental  Observations 

Slip  lines  are  often  developed  in  the  early  stage  of  fatigue  cycling  to  produce  ’'soft" 
regions  where  local  deformation  tends  to  concentrate  as  shown  by  Woods  (1973),  Winter 
(1974),  Finney  and  Laird  (1975).  Slip  lines  generally  multiply  and  grow  into  persistent  slip 
bands.  Persistent  slip  bands  carry  essentially  all  the  plastic  strain  in  low  amplitude  fatigue 
tests  as  indicated  by  Brown  and  Ogin  (1984).  Tests  on  aluminum  single  crystals  under  cyclic 
tension  and  compression  by  Charsley  and  Thompson  (1963)  have  shown  that  a  reversal  of 
stress  after  a  prior  tensile  deformation  gives  rise  to  new  slip  lines.  Similar  "compressive"  slip 
lines  were  observed  by  Buckley  and  Entwistle  (1956)  to  form  between  "tensile"  slip  lines. 
Forsyth  (1954)  found  that  the  slip  lines  produced  in  the  forward  loading  and  those  produced 
in  the  reversed  loading  were  in  the  same  fatigue  band,  very  close  but  distinct  from  each  other. 
Wood  and  Bendler  (1962)  have  shown  with  their  electron-micrographs,  that  both  forward  and 
backward  slips  have  occurred  within  a  fatigue  band;  also  the  slip  that  occurs  in  cyclic  loading 
does  not  lead  to  significant  deformation  in  the  bulk  of  the  matrix.  This  is  verified  by  the  lack 
of  dark  spots  on  the  X-ray  reflection  photographs  of  fatigue  specimens  shown  by  Woods 
(1956).  Intrusions  are  often  considered  as  embryos  of  cracks.  Other  experimental  observa¬ 
tions  and  suggested  models  of  fatigue  crack  initiation  were  given  in  two  excellent  reviews  by 
Grosskreutz  (1971)  and  Laird  and  Duquette  (1972).  These  experimental  observations  show 
that  the  slip  lines  formed  in  the  forward  loading  and  those  formed  in  the  reversed  loading  in 
fatigue  specimens  often  are  closely  spaced  but  distinct  from  each  other.  Observations  clearly 
support  the  proposed  gating  mechanism  provided  by  the  microstress  fields  due  to  two  alterna¬ 
tively  sliding  slices. 
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Extent  of  Extrusion  and  Intrusion: 


The  buildup  of  the  slip  strain  ea}  in  P  and  Q  is  caused  by  e„a  in  R.  If  R  were  cut 
out,  the  free  length  of  R  would  be  longer  than  the  slot  by  an  amount  referred  to  as  the  "static 
extrusion"  by  Mughrabi  et  al.  (1983).  This  e^  causes  an  initial  compression  x^  in  R. 
Under  cyclic  loading,  the  extrusion  grows  and  the  thin  slice  R  increases  in  length.  This  elon¬ 
gation  causes  the  compression  to  decrease.  A  question  raised  by  Mughrabi  (1980)  and 
Essmann,  Gosele  and  Mughrabi  (1981)  is  whether  the  extrusion  growth  will  cease  after  the 
extrusion  has  reached  the  magnitude  of  the  static  extrusion.  This  question  is  to  be  discussed 
in  the  following  calculations. 

There  are  twelve  slip  systems  in  a  f.c.c.  crystal.  The  change  of  the  direct  stress  xaa  in 
R  causes  changes  of  resolved  shear  stresses  in  all  slip  systems.  When  the  decrease  of 
compression  R  becomes  large,  its  resulting  residual  stresses  combined  with  the  applied  stress 
can  cause  a  second  slip  system  to  have  shear  stress  reaching  the  critical  and  slide.  The  plastic 
strain  e^,  caused  by  slip  in  this  slip  system  has  a  tensor  component  eaa  just  like  e^  in  caus¬ 
ing  the  positive  and  negative  x^p  in  P  and  Q  respectively.  Hence  this  secondary  slip  can 
increase  greatly  the  extent  of  extrusion  and  intrusion. 

Method  for  Calculating  Incremental  Slip  in  P,  Q  and  R: 

The  slip  direction  a  and  the  normal  to  the  slip  plane  [3  make  45°  to  the  loading  axis 
(Fig.  1).  As  the  extrusion  grows,  ea„  in  the  slip  band  causes  xaa  to  increase  and  can  cause  a 
second  slip  system  to  slide.  aP  lie  in  the  x1?  x2  plane  as  shown  in  Fig.  1  but  ^r|  is  not 
confined  to  the  xl5  x2-plane.  Hence  the  equivalent  force  F^  due  to  es-  can  have  component 
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normal  to  this  plane.  In  1972,  Lin  and  Lin  showed  the  effect  of  a  secondary  slip  on  fatigue 
crack  initiation.  To  simplify  the  numerical  computation,  they  assumed  the  second  slip  system 
to  be  on  the  xrplane  along  the  x2  direction  or  vice-versa.  This  is  not  an  actual  slip  system. 


In  the  present  study,  the  second  slip  system  analyzed  is  one  actually  existing  in  the  crystal. 
The  presence  of  a  body  force  component  along  the  x3  direction  requires  the  modification  of 
the  solution  of  the  plane  strain  problem.  A  similar  problem  was  shown  in  the  analysis  of 
prismatic  anisotropic  bars  by  Lekhnitski  (1963)  and  is  referred  to  as  the  generalized  plane 
strain  problem.  This  generalized  plane  strain  is  expressed  as 


Uj  =  Uj(x j  x2),  i  =  1,2,3 


This  gives 


*ij  =  2G 


5,j  0  +  >/2  (u^  +  uj  j) 


where  9  =  uJ(1  +  u2>2  and  v  is  Poisson’s  ratio.  The  substitution  of  this  expression  into  the 
condition  of  equilibrium  (3)  yields 


V2  U“  +  7^2v  +  "g"  =  °  ’  a=U 


V2  u,  +  -f  =  0 


The  differential  equations  (9)  and  (10)  are  not  coupLd  and  can  be  solved  separately. 


The  three  slip  directions  of  each  of  the  four  slip  planes  of  a  f.c.c.  crystal  are  shown  in 
Fig.  2.  The  most  favorably  oriented  crystal  of  a  polycrystal  loaded  under  alternate  tension 
and  compression  along  the  x2  axis  (Fig.  1)  has  a  slip  plane  and  a  slip  direction  making  45° 
with  the  direction  of  loading.  The  slip  system  is  referred  to  as  the  primary  slip  system.  Let 
a-,  axis  (Fig.  2),  correspond  to  this  system.  During  fatigue  loading,  the  build-up  of  large  local 
plastic  shear  strain  in  the  primary  slip  system,  positive  in  P  and  negative  in  Q,  tends  to  start 
an  extrusion  or  an  intrusion.  Consequently  an  appreciable  direct  stress  x ^  will  be  built-up  in 
R.  The  Schmid’s  factors  of  all  the  12  slip  systems  for  the  stress  x,j  and  those  for  xaot  are 
listed  in  Table  1.  It  is  seen  that  there  are  four  slip  systems  Cj,  c3,  dj  and  d2  equally  favorable 
under  xaa.  Among  these  four,  c3  is  most  favorably  oriented  under  x22  where  x22  is  the  alter¬ 
nate  loading.  Hence  c3  is  considered  to  be  the  active  second  slip  system  in  R. 


Plastic  strain  eap  due  to  slip  in  the  primary  slip  system  ccp,  causes  equivalent  forces  Fq 
and  F2.  This  is  a  plane  strain  problem  Melan  has  given  the  solution  of  the  stress  distribution 
caused  by  a  point  force  in  a  semi-infinite  plate  in  plane  stress.  Melan’s  solution  was  modified 
for  plane  strain  by  Tung  and  Lin  (1966).  Let  x,j  (x,  x)  be  the  stress  at  x  due  to  a  unit  force 
along  xk  axis  applied  at  x  and  q>k  (x,  IT)  be  the  corresponding  stress  function,  where  x  denotes 
(xj,  x2,  x3).  Lin  and  Lin,  1974  has  expressed  the  stress  components  in  terms  of  the  stress 
functions  as: 


k  ,  _  32*k 

til  (x,  x)  =  TT  ■ 

C/Xt“ 


.  326k 

X22  (X,  X)  =  , 

ox  i* 


d2<i>k 

tf*  <*■  *>  =  '  k=u 


(H) 


where  the  stress  functions  are 


i 
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Slip  Planes 


Slip  Directions 


Plane  ^  =  -|  [  1  1  0  ]  c,  =  ~  [  1  1  0  ] 

ABE  (111)  a2  =  -—  [  T  0  1  ]  c2  =  —  [  1  0  1  ] 

2  2 

BCE  (ill)  a3  =  y  [  0  1  T  ]  c3  =  —  foil) 

2  1 

CDE(TTl)  bj  =  -J  [  T T 0  ]  d,  =  -—  [  1  1  0  1 

DAE  (Til)  b2  =  ~  [  0  1  1  ]  d2  =  ~  t  0  T  1  ] 

b3  =  J  [  1  0  T  ]  d3  =  |  [  7  0  T  1 

Agure  2.  Crystallographic  Directions  of  a  f.c.c.  Crystal 


Stresses  in  Different  Slip  System  Caused  by  Cyclic 
C22  and  xlut  (a  denotes  a2  Direction) 


-0.197  -0.272  -0  469  0.469  -0.197  0.272 


Plastic  strain  due  to  slip  in  the  secondary  slip  system  causes  F3  in  addition  to  F] 
and  Fi  To  solve  Eq.  (10)  for  F3,  we  write  <>3  (x,  x)  as  G  u3(x,  x).  The  stress  is  then 

3  ,  _  3  ,  __ 

xf3  (x,  X)  =  —  ;  t23  (X,  X)  =  — 

dxi  dx2  (i4) 

All  other  stress  component’s  are  zero.  For  a  unit  concentrated  force  F3  at  x  Eq.  10  gives 
V2  <{>3  +  5  (x  -  x)  =  0 

where  5(x,  x)  is  the  Dirac  Delta  function.  With  the  boundary  conditions  of 
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9<}>3 

— —  =  0  at  x,  =0  and 

o'-, 


a<^  ,  303  r~  7 

- —  and  — —  — >  0  as  \xi  +  x-r  — >  , 

oxj  oxi 


we  have 


03(x,  x)  =  -|s(/n  X^/nXj) 


where  X[  and  X2  are  given  in  Eq.  13. 


(15) 


(16) 


Numerical  Calculation  and  Results: 

The  thicknesses  of  the  slices  P,  Q  and  R  are  very  small  as  compared  to  their  lengths. 
As  discussed  before,  the  slip  strains  in  these  slices  are  taken  to  be  constant  across  the  thick¬ 
ness.  The  slices  are  divided  into  a  number  of  parallegram  grids  along  the  length  in  the  a 
direction.  For  numerical  calculation,  a  grid  of  constant  plastic  strain  has  been  used.  The 
corresponding  grid  stress  is  the  average  stress  in  the  grid.  From  the  plane  strain  solution  of  a 
semi-infinite  medium  the  stress  field  caused  by  a  uniform  plastic  strain  eaPn  in  the  n-th  grid 
was  calculated.  The  relief  of  shear  stress  t£p  (x)  at  any  point  is  expressed  as 


(x)  =  -  C  (x,  ap  ;  n,  ap)  e”Pn 


> 

(17) 


where  C(x,  aP  ;  n,  ap)  is  the  resolved  shear  stress  at  x  due  to  unit  uniform  plastic  strain  eap 
in  grid  "n".  The  average  relieved  stress  x„p(x)  over  the  m-th  grid  denoted  by  taPm  and  is 


written  as 


I 
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=  -  C  (m,  ap  ;  n  ,  ap)  eaPc 


As  a  cyclic  loading  proceeds,  slips  occur  in  P  and  Q  and  cause  a  residual  stress  field. 
The  total  resolved  shear  stress  in  the  m111  gnd  of  the  a(3  slip  system  is  the  sum  of  the  initial, 
residual  and  applied  stresses: 


=  Tip.  +  -  X  C  (m,  ap  ;  n,  ap)  e^ 

n 

where  "n"  is  summed  over  all  the  sliding  grids  with  plastic  strain  eapn. 


(18) 


Plastic  deformation  in  metals  is  highly  localized  (Brown,  1952).  Plastic  strain  concen¬ 
trates  in  thin  slices.  The  microscopic  plastic  strain  in  the  thin  slice  is  much  larger  than  the 
macroscopic  plastic  strain  in  the  metal.  Hence,  the  strain-hardening  of  the  thin  slice  is  much 
less  than  that  of  the  metal  and  is  neglected.  This  gives  a  constant  critical  resolved  shear 
stress  xc.  Sliding  occurs  in  the  grids  where  T  (x)  =  ±  Tc.  For  an  incremental  applied  shear 
stress  A  x£,  Axm  =  At  =  0.  Eq.  1 8  yields 


Kp.  =  X  c  ap  ;  n-  ap)  Aeaj}n 

n 


(19) 


There  are  as  many  unknown  Aeapn’s  as  the  number  of  the  above  equations.  The  plastic  strain 
increments  Aeapn  in  the  sliding  grids  for  an  incremental  loading  Ata  can  be  readily  calculated. 
Similarly,  the  stress  components  x„a  and  Xpp  in  R  can  be  calculated  (Lin,  1972)  and  written  as 

=  X  c  (m-  aa  ;  n,  ap)  ea'Pn 


(20) 
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if  pm  =  L  c  (m,  PP  ;  n,  ctp)  e  "pn 

n 

From  these  calculations,  it  was  found  that  Xppm  in  R  is  quite  small. 


To  find  C  (m,  aa  ;  n,  aP),  a  different  approach,  which  gives  a  better  physical  insight, 
is  considered.  The  elastic  shear  strain  eap,  that  equals  xc/2G  in  the  sliding  grids,  is  much  less 
than  the  cumulative  plastic  strain  eap.  The  cumulative  plastic  shear  strain  eap  in  P  and  Q  are 
about  the  same.  Denoting  the  thickness  of  the  slices  P  and  O  by  t,  we  can  express  the  dis¬ 
placement  in  R  along  the  a  direction  as  -2teaj3pg  or  approximately  -2teappg.  The  negative 
sign  is  due  to  the  fact  that  the  direction  of  the  extrusion  is  opposite  to  the  direction  of  the  a 
direction.  This  causes  a  tensile  strain  in  R, 


3e„ 


e  =  -  2t  — ^ 
e“aR  1  3a 


(21) 


Since  Xpp  is  very  small  and  assumed  to  be  zero,  we  have  the  residual  tensile  stress 


Tr  =: 
'rtrt 


3c, 


l_v2  da 


(22) 


With  the  residual  stress  x^  increasing  with  cycles  of  loading,  the  initial  compressive  stress 
laa  is  gradually  relieved. 


As  discussed  earlier,  C3  (Table  I)  is  the  most  likely  second  slip  system  to  become 
active.  The  resolved  shear  stress  in  c3  is  denoted  by  Xe^.  As  shown  in  Table  1,  the  initial 
resolved  shear  stress  x^  varies  with  the  initial  stress  components  xj 3,  x33,  x33.  For  our 
numerical  calculation  x^  is  assumed  to  be  zero.  From  Table  1, 
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A  x^  =  -  .318  A  x22  -  0.408  A  x^ 

When  x^  increase  to  x°  or  decreases  to  -xc,  this  slip  system  slides  and  causes  et. 
have  slip  in  the  ap  slip  system  in  P  and  Q  and  in  the  ^T)  slip  system  in  R. 


(23) 
Now  we 


Let  the  two  slices  P  and  Q  be  divided  into  2N  grids  and  R  into  M  grids.  The  resolved 
shear  stress  in  the  ap  slip  system  of  the  i-th  grid  is  written  as 

x£p.  =  -  C(i,  ap  ;  j,  ap)  e"Pj  -  C(i,  ap  ;  k,  ^ri)  e^ 

x^_  =  -  C(i,  ^  ;  j,  ap)  e"Pj  -  C(i,  ;  k<  &\)  e^k 

The  repetition  of  subscript  "j"  denotes  summation  from  one  to  2N  and  that  of  subscript  "k" 
denotes  summation  from  one  to  M.  C(i,  £r|  ;  j,  ap)  is  the  residual  stress  x^  in  i-th  grid  due 
to  a  unit  plastic  strain  eap  in  the  j-th  grid.  For  constant  strain  erap  in  grid,  j  the  corresponding 
C(i,  ^r|  ;  j,  aP)  is  the  average  shear  stress  x^  in  the  i-th  grid.  A  method  of  numerical  calcu¬ 
lation  using  grids  with  strain  linear  in  a  direction  has  also  been  developed  and  used.  With 
residual  stresses  given  by  Eq.  (24),  the  incremental  plastic  strain  in  all  grids  are  calculated  by 
Eq.  (19). 


The  dimensions  of  the  most  favorably-oriented  crystal  at  the  free  surface  are  shown  in 

Fig.  1.  The  thickness  in  P  and  Q  is  0.01  p  along  x2  dimension  and  thickness  in  R  is  0.1  p 

along  the  same  direction.  The  initial  strain  c^a  in  R  is  assumed  to  vary  linearly  from  zero  at 

the  free  surface  to  3.64  x  10-5  at  Xj  =  50|i.  The  static  extrusion  is  then  Vz  x.  3.64  x  10~5  x 

r-  ,  dua  „  . 

50  x  v2  p  =  1.28  x  10  p.  ea«  =  Vz  .  This  corresponds  to  a  plastic  strain  ean  in  P  and 

H  dp 

Q  at  the  free  surface  of  0.091. 
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This  initial  strain  field  causes  a  uniform  initial  shear  stress  of  1.46  x  1CT3  MPa  (0.20 
p.s.i.)  in  P  and  -1.46  x  10-3  MPa  (0.20  p.s.i.)  in  Q.  The  polycrystal  is  subject  to  a  cyclic  ten¬ 
sile  and  compressive  loading  in  the  x2  direction.  Let  xE  be  the  excessive  shear  stress  defined 
as  the  initial  shear  stress  x1  plus  the  maximum  applied  shear  stress  xa  minus  the  critical  shear 
stress  xc: 

xE  =  x1  +  xa  -  xc 

(25) 

Tests  on  single  aluminum  crystals  under  cyclic  loading  under  constant  plastic  strain  amplitude 
by  Snowden,  1963,  show  that  the  resolved  shear  stress  increases  in  the  first  50  cycles  and 
approaches  a  saturated  value  of  about  2.94  MPa  (425.8  p.s.i.  or  300gm/nm2).  Hence  xc  is 
taken  as  2.94  MPa  (425.8  p.s.i.).  For  a  maximum  applied  shear  stress  xa  of  2.94  MPa  (425.8 
p.s.i.),  the  excessive  shear  stress  xE  is  1.46  x  10-3  MPa  (0.20  p.s.i.). 

With  the  process  of  sliding  discussed  earlier,  the  plastic  strain  distributions  e^  in  the 
secondary  slip  system  are  shown  in  Fig.  3  and  the  corresponding  resolved  shear  stress  (x^) 
variations  are  shown  in  Fig.  4.  The  plastic  strain  in  the  second  slip  system  causes  some 
eaa<  which  in  turn,  causes  positive  Xap  in  P  and  a  negative  xap  in  Q  just  like  the  initial 
Strain  e ^ .  This  e^  causes  additional  +xap  in  P  and  -xap  in  Q.  The  sum  of  the  initial  plus 
residual  shear  stresses  x^p  +  x£p,  in  P  and  Q,  on  which  the  extrusion  growth  depends,  is 
shown  in  Fig.  5.  It  is  seen  that  this  sum  is  considerably  more  in  the  cases  with  the  secondary 
slip  than  in  those  without.  The  plastic  shear  strain  distributions  eap  in  P  and  Q  at  different 
cycles  of  loading  were  calculated  for  the  cases  with  and  without  the  secondary  slip  and  are 
shown  respectively  as  solid  and  dotted  lines  in  Fig.  6.  The  ejp’s  at  the  free  surfaces  vs. 


2500  Cycles 


ic  Strain  Distribution  of  the  Second 
Slip  System 
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Residual  Stress  Distribution  in  the  Se 
Slip  System 


cycles  of  loading  are  shown  in  Fig.  7.  It  is  seen  for  the  case  without  considering  the  secon¬ 
dary  slip,  extrusion  growth  ceases  as  the  extrusion  approaches  the  static  extrusion  The  plastic 
strains  eap  in  P  and  Q  at  the  free  surface,  which  represent  the  amount  of  extrusion  or  intru¬ 
sion,  with  the  secondary  slip  is  four  times  that  without  the  secondary  slip. 

Conclusions 

Under  cyclic  tension  and  compression  the  plastic  shear  strain  growth  eap  in  P  and  Q 
were  calculated  with  and  without  considering  the  secondary  slip  system  ^r[.  The  plastic  strain 
c'Sti  in  the  secondary  slip  system  has  a  component  eaa  which  has  the  same  effect  in  causing 
xa|}  in  ?  and  Q  as  e‘aa  ■  Hence  ean  in  P  and  Q  with  the  secondary  slip  in  the  present  case  is 
four  times  more  than  that  without  the  secondary  slip.  Consequently,  the  extent  of  extrusion  or 
intrusion  in  the  cases  with  the  secondary  slip  is  greatly  more  than  those  without.  This 
explains  the  observed  extrusion  growth  beyond  the  static  extrusion  as  reported  by  Mughrabi 
(1983)  and  explains  that  a  face-centered  polycrystal  has  greater  extrusion  height  and  intrusion 
depth  than  hexagonal  polycrystals.  The  stress  intensity  factor  of  crack  depends  on  the  extent 
of  extrusion  and/or  intrusion.  Hence  this  extent  is  important  for  the  study  of  crack  initiation. 
The  initial  strain  distribution  e„a  assumed  in  the  slice  R  is  arbitrary.  Similar  calculations  can 
be  made  for  other  initial  strain  distributions.  The  present  method  of  analysis  satisfies  the 
equilibrium,  continuity  of  displacement  and  slip  characteristics  of  the  crystals. 


50  ( 


(~)D  pirn  (+)d  J> 
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itial  and  Residua)  Stress  Distribution 
’  and  Q  (P  and  Q  opposite  sign) 


Acknowledgments 


The  support  of  the  U.S.  Office  of  Naval  Research  through  Contract  N00014-86-K-0153 
and  the  interest  of  the  scientific  officer.  Dr.  Yapa  Rajapakse.  are  gratefully  acknowledged. 


26 


References 


1.  Antonopoulus,  J.G.,  Brown,  L.M.,  and  Winter,  A.J.,  1976,  Phil.  Mag.  34,  549-563. 

2.  Brown,  A.F.,  1952,  Adv.  Phys.  1,  427. 

3.  Brown,  L.M.  and  Ogin,  S.L.  1984,  "Fundamentals  of  Deformation  and  Fracture," 
Eshelby  Memorial  Symposium  500-528. 

4.  Buckley,  S.N.  and  Entwistle,  K.M.,  1956,  ACTA  Metall.  4,  352-361. 

5.  Charsely,  P.  and  Thompson  N,  1963,  Phil.  Mag.  8,  77-85. 

6.  Cottrell,  A.H.  and  Hull,  D.  1957,  Proc.  Rog.  Soc.  A.  202-211. 

7.  Eshelby,  J.D.  1957.  Proc.  Roy.  Sec.  A  241  pp.  376-396. 

8.  Essmann,  V.,  Gosele  V.  and  Mughrabi  H.  1981,  Phil,  Mag,  44.  405-426. 

9.  Finney,  J.M.  and  Laird  C.  1975,  Phil.  Mag.  31,  339-366. 

10.  Forsyth,  P.J.E.  1953  Nature  171,  172. 

11.  Forsyth,  P.J.E.  1954,  J.  Inst.  Metals  82,  449-454. 

12.  Grosskreutz,  J.C.  1971  Phys.  Stat.  Solids  (b)  47,  359-396. 

13.  Hull,  D.  1958,  J.  Inst,  Metals,  86,  425. 

14.  Kennedy,  A.J.  1963  "Process  of  Creep  and  Fatigue  of  Metals",  John  Wiley  &  Sons, 
pp.  331-346. 

15.  Laird,  E.  C.  and  Duquette  D.J.  1972  "Corrosion  Fatigue"  Edited  by  Devereux  O.F., 
McEviley,  A.J.  and  Stachle,  R.W.  pp.  88-117 

16.  Lin,  T.H.  and  Ito,  Y.M.  1967  J.  Appl.  Phys.  3,  775-780. 

17.  Lin.  T.H.  1968  "Theory  of  Inelastic  Structures"  John  Wiley  &  Sons.,  pp.  48-49. 

18.  Lin,  T.H.  and  Ito,  Y.M.  1969  J.  Mech.  Phys.  Solids  17,  511-527. 

19.  Lin,  S.R  1972  Ph.D.  Thesis,  Univ.  Calif.  Los  Angeles. 

20.  Lin,  S.R.  and  Lin,  T.H.  1974  J.  Mech.  Phys.  Solids,  22,  177-192. 

21.  Lin,  T.H.  1977  Rev.  Deformation  Behavior  Materials,  Fdited  by  P.  Feltb'im.  Freund 

Publishing  House,  Israel,  264-314. 


27 


22.  Lin,  T.H.  and  Lin,  S.R.  1983.  J.  Appl.  Mech.  50,  367-372. 

23.  McEviley  J.  Jr.  and  Machlin  E.S.  1959  Proc.  Int.  Conf.  Atomic  Mechanism  of  Frac¬ 
tion  Technology,  John  Wiley  &  Sons. 

24.  Meke,  K.  and  Blochwitz,  C.  1980  Phys.  Stat.  Sol.  a  61,  5. 

25.  Mott,  N.F.  1958  ACTA  Metall,  6,  195. 

26.  Mughrabi,  H.  1980  Microscopic  Mechanisms  of  Fatigue  Strength  of  Metals  and 
Alloys.  Vol.  3  Perganon  Press. 

27.  Mughrabi  H.,  Wang  R.,  Differt  K.  and  Essmann  V.  1983,  Fatigue  Mechanism,  ASTM 
STP  811,  p.  5-45. 

28.  Parker,  E.R.,  1961  "Mechanical  Behavior  of  Materials  at  Elevated  Temperatures” 
Edited  by  Dorn  J .£.,  McGraw-Hill  Book  Co.  pp.  129-148. 

29.  Snowden,  K.V.,  1963,  ACTA  Metallugica  Vol.  11  No.7,  675-684. 

30.  Thompson,  N.  Wadsworth  N.J.  and  Louat  N.  1956,  Phil.  Mag.  1,  113. 

31.  Thompson,  N.  1959  Proc.  Int.  Conf.  Atomic  Mechanism  of  Fraction  Technology,  John 

Wiley  &  Sons. 

32.  Winter,  A.T.  1974,  Phil.  Mag.  30,  719-738. 

33.  Wood,  W.A.  1956  "Fatigue  in  Aircraft  Structures",  Edited  by  A.F.  Freudenthal. 
Academic  Press,  New  York.  pp.  1-19. 

34.  Wood,  W.A.  and  Bendler,  A.M.  1962,  Trans.  Metall.  Soc.  AIME  224  180-186. 

35.  Wood,  P.J.  1973,  Phil.  Mag.  30,  719-738. 


28 


